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Abstract. Andredakis, Peletier & Balcells (1995) fit Ser- 
sic's law fi(r) ~ r 1 /" to the bulges of the Balcells & 
Peletier (1994) galaxy sample, and infer that n drops with 
morphological type T from n rs 4-6 for SO to n = 1 (expo- 
nential) for Sc's. We use collisionless N body simulations 
to test the assumption that initially the surface bright- 
ness profiles of all bulges were exponential, and that the 
steepening of the profiles toward the early-types is due to 
satellite accretion. The results are positive. After the ac- 
cretion of a satellite, bulge-disk fits show that the bulge 
grows and that the bulge profile index n increases propor- 
tional to the satellite mass. For a satellite as massive as 
the bulge, n rises from 1 to 4. We present kinematic di- 
agnostics on the remnants and disk thickening. The latter 
suggests that the bulge growth must have occurred before 
the last formation of a thin disk in the galaxy. The thick 
disks created by the merger are reminiscent of thick disks 
seen in early-type edge-on galaxies. The efficiency of the 
process suggests that present day bulges of late-type spi- 
rals showing exponential profiles cannot have grown sig- 
nificantly by collisionless mergers. 
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1. Introduction 

Two paradigms are commonly used to explain the forma- 
tion of the central bulges of disk galaxies. In the first, the 
bulge forms prior to the disk, perhaps as part or an after- 
math of the formation of the stellar halo (Gilmore & Wyse 
1998), or as a result of early merging resulting in an r 1 / 4 
spheroid around which a new thin disk forms from sur- 
rounding gas (eg. Kauffmann, Guiderdoni & White 1994). 
In the second, bulges form and grow after bar instabili- 
ties draw disk stars and gas to the center (Pfenniger & 
Norman 1990). Details about these paradigms are given 
in the review by Wyse, Gilmore & Franx (1997). Bouwens, 
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Cayon & Silk (1999) perform a comparative test of these 
hypotheses. 

What is the role of interactions, mergers and accretion 
of satellites in the functioning of these processes?. Satel- 
lite accretion must have occurred several times in a disk 
galaxy over a Hubble time. While mergers of galaxies with 
similar mass fully destroy the disk (eg. Barnes & Hernquist 
1991) and leave remnants resembling an elliptical galaxy, 
satellite accretion does less damage to the disk; and, if 
the satellite is dense enough, it must reach the galaxy 
center. The accretion in this case may drive the growth of 
the central bulge. Several studies have addressed whether 
mergers are responsible for pushing a spiral galaxy along 
the Hubble sequence toward Sa/SO (eg. Walker et al 1996). 
Mergers must have been efficient in groups and in the field; 
in the cluster environment, disk cosmological fading due 
to the evaporation of neutral gas, the cluster tidal field 
(Dubinski 1998) and occasional galaxy-galaxy encounters 
(Moore et al. 1996) contribute as well to the growing 
prominence of the bulge. 

Satellite accretion onto disk galaxies has been exten- 
sively studied in the past using N body techniques. The 
focus of these studies have been the thickening of the disk 
(Quinn & Goodman 1986, Toth & Ostriker 1992, Quinn, 
Hernquist & Fullagar 1993, hereafter QHF93, Walker et al 

1996, Huang & Carlberg 1997), the formation of counter- 
rotating disks in spirals (Thakar & Ryden 1996), the for- 
mation of X structures in SO galaxies (Mihos et al. 1995), 
and the triggering of nuclear starbursts (Mihos & Hern- 
quist 1994). Diagnostics related to the central bulge are 
largely ignored. Indeed, except for the latter paper, none 
of the studies mentioned above include a central bulge in 
their initial galaxy model, a bold simplification given that 
the potential of the central bulge is a key ingredient of the 
merger dynamics of both stars (Hernquist 1993) and gas 
(Mihos & Hernquist 1994). 

Sufficient data exists nowadays on bulges (Balcells & 
Peletier 1994, hereafter BP94, Peletier & Balcells 1996, 

1997, Peletier et al. 1999, Carollo 1999, de Jong 1996, 
Prieto et al. 1999, Ratnatunga, Griffiths & Ostrander 

1998, Schade et al. 1996, Abraham et al. 1999, Marleau 
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& Simard 1998) to allow the structure and dynamics of 
bulges to be used as constraints on the accretion history 
of disk galaxies. However, no detailed studies of the growth 
of bulges by accretion of satellites have been published to 
date. 

In this paper we study whether the shape of the sur- 
face brightness profiles of bulges changes with the merger 
history of spiral galaxies. Andredakis, Peletier & Balcells 
(1995, hereafter APB95) show that bulges follow a system- 
atic behavior in the n vs. T plane, where n is the exponent 
index in Sersic's law \x(r) ~ r 1 '™, and T is the galaxy type 
index given eg. in the RC3. APB95 find that n is 4-6 in 
early, SO types, decreasing to n » 1 for late types (Sc). 

The APB95 result could result from the effect of the 
disk potential on a bulge initially following an r 1 / 4 law, 
in an scenario in which the bulge forms prior to the disk. 
Andredakis (1998) studies the adiabatic growth of a disk 
onto an existing r 1 / 4 spheroid, and finds that the disk 
potential does modify the bulge surface brightness profile, 
lowering the n in the exponent of Sersic's law. But the 
mechanism saturates at n = 2, showing that exponential 
bulges (n = 1) cannot be produced by adiabatic growth 
of the disk around an r 1 / 4 bulge. 

That merging is an efficient way to generate r 1 / 4 
surface brightness profiles is well-known (Gerhard 1981; 
Barnes 1988). We may thus conjecture that the behavior 
found by APB95 reflects the effects of the relaxation in 
an scenario of bulge growth by accretion. Here we ad- 
dress the following question: assuming that bulges are 
formed with exponential profiles, does the accretion of 
dense satellites simultaneously drive the growth of the 
bulge mass and the evolution of the bulge surface density 
profile toward higher-n Sersic law shapes? In cosmological 
models of galaxy formation involving gas and stars, cen- 
tral bulges show characteristic exponential surface density 
profiles (Domfnguez-Tenreiro et al. 1999), suggesting that 
bulges may start out with exponential profiles, and that 
any changes are imprints of subsequent evolution. To test 
the bulge growth conjecture we consider disk-bulge-halo 
galaxy models for which the bulge surface brightness pro- 
file is initially exponential, and fit Sersic's law to the bulge 
surface brightness profile after an accretion of a satellite. 
We then draw "growth vectors" in the n vs. T plane. 

Our results suggest that accretion of dense satellites is 
efficient in turning exponential bulges into r x / n>1 bulges. 
After the merger, the n index of Sersic's law increases 
proportional to the satellite mass. It reaches n — 4 for a 
satellite as massive as the original bulge. Our models sug- 
gest that this result is largely independent of the internal 
structure of the satellite, as long as the satellite is dense 
enough to reach the center undisrupted. We quantity disk 
thickening. We concur with other studies of spiral galaxy 
evolution via merging, that only if the thin disk is rebuilt 
after the accretion can we match real galaxies. We analyze 
the kinematic structure of the merger remnants to provide 
additional tests of the model. For direct mergers, the ro- 



tation curve shapes are steeper for more massive bulges, 
consistent with observations. Massive satellites on retro- 
grade orbits result in counterrotating bulges, suggesting 
either that massive accretion events are rare, or that ac- 
cretion preferentially occurs on prograde orbits. 

Model details are given in § 0. The results on the shape 
of the bulge surface brightness profiles are given in § |[ 
Model kinematics are presented in § [|. Section || briefly 
presents results on disk heating. Implications are discussed 
in § ^[ A summary of results is given in § ^. Throughout 
the paper, M denotes mass and not absolute magnitude. 

2. Models 

The primary galaxy is modeled with the bulge-disk-halo 
model of Kuijken & Dubinski (1995, hereafter KD95). The 
bulge is modeled as a King model. This model was chosen 
for the bulge because its surface density profile is expo- 
nential between the core radius and the truncation radius, 
providing a good match to the surface brightness profiles 
of late-type spirals, and because it is spheroidal. The core 
radius is 0.15 (units given below), and the concentration 
parameter, taken as the ratio of the bulge truncation ra- 
dius and the core radius, is 6.7. The disk surface density is 
exponential both in the galactic plane and in the perpen- 
dicular direction to the plane. The scale height is ten times 
less than the scale length (see eg. Guthrie 1992, de Grijs 
1998). The chosen disk velocity dispersion makes the disk 
warm with a Toomre Q — 1.7 at the disk half mass ra- 
dius. The value of Q is fairly constant throughout the disk, 
though rising both in the center within 1 scale length and 
near the edge (KD95). This Q parameter does not allow 
the growing of bar-type perturbations in the disk when in 
isolation. The halo has a distribution function of an Evans 
model (Kuijken & Dubinski 1994). 

The model used for our experiment is equal to model 
A of KD95 This model matches the Milky Way when the 
units of length, velocity and mass are R = 4.5 kpc, V — 
220 km s -1 , M = 5.1 x 10 10 Af Q . Note that the core radius 
is then 0.71 kpc, indicating that the models lack central 
resolution. Masses, radii, and number of particles for each 
component are given in Table [l]. A gravitational constant 
of G — 1 is used throughout. 

The satellite galaxy is modeled either as a non-rotating 
Hernquist sphere (Hernquist 1990), with outer radius sim- 
ilar to that of the bulge, or with a King model. The King 
model satellite is useful to verify that the accretion-driven 
evolution toward an r 1 / 4 profile is not motivated by the 
fact that the Hernquist satellite models already have a 
surface density profile approaching the r 1 / 4 law. 

Both the disk and the satellite galaxy models are al- 
lowed to relax separately for 1.5 disk rotation periods 
(measured at 1 disk scale length) before starting the sim- 
ulations. The surface density profile of the initial disk 
galaxy model after relaxation is shown in figure [|a. Its 
projected rotation curve is shown in figure [7ja together 
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Table 1. Initial parameters of the main galaxy 
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Fig. 1. (See aguerri.figl.gif file) The time evolution of the luminous matter for model 1. Views from the z, x and y 
axes are given. At t = 0, the satellite is at (x,y,z) = (8.66, 0., 5.). Orbit is counterclockwise, prograde with the disk 
rotation. 



Fig. 2. (See aguerri.fig2.gif file) The time evolution of the luminous matter for model 5. Views from the z, x and y 
axes are given. At t = 0, the satellite is at (x,y,z) = (8.66, 0., 5.). Orbit is counterclockwise, prograde with the disk 
rotation. 



with the model's circular velocity curve. We ran the disk 
galaxy model in isolation for 2000 time updates (14 disk 
rotation periods at 1 disk scale length, the maximum du- 
ration of our merger experiments) . The shape of the bulge 
surface density profile is stable (same n value at start and 
at end of run), while the disk thickens in close agreement 
with the thickening results of KD95. This demonstrates 
that, while the number of particles is admittedly low, the 
results do not suffer from two-body relaxation as the more 
massive halo particles cross the disk and the bulge. This 
result also demonstrates the excellent stability properties 
of the KD95 model, which makes it ideally suited to ex- 
periments of interaction-driven structural evolution such 
as the ones in this paper. 

We run several merger experiments varying the mass 
of the satellite galaxy. We explored three values of the 
bulge-to-satellite mass ratio, 1:1, 3:1, and 6:1. For each 
satellite mass setting we run a direct orbit (inclination 30° 
w.r.t. the disk plane) and a retrograde orbit (inclination 
150°). Initial orbits were elliptical with apocenter equal to 
twice the disk outer radius and pericenter equal to twice 
the diskscale length. We don't explore the dependency of 
the results on the orbital energy or angular momentum. 
This choice is probably restrictive. We expect a depen- 
dency on orbital initial conditions to be significant for the 
disk (eg. QHF93), though possibly not for the bulge; when 
the satellite reaches the inner parts of the galaxy, the po- 
tential's circular velocity should determine the energy of 
the bulge-satellite merger more so than the orbit initial 
conditions. For all these models a Hernquist satellite was 
used. Retrograde merger orbits with mass ratios 1:1 and 
1:3 were repeated using a King model satellite. While full 
exploration of satellite density effects is beyond the scope 
of this paper, we run one model in which the satellite has 



low density (model 7 in Table |J) . This model is similar to 
the models of Velazquez & White (1999). Orbital param- 
eters for the merger experiments are shown in Table 2. 
Satellite masses and half-mass radii are given in Table 2. 
We used a M ~ r 1 ' 3 scaling between satellites of different 
masses. In the remainder of the paper, models are referred 
to using a three-character code: SM0, where S describes the 
satellite (H for Hernquist, K for King, L for low-density), 
M describes the ratio of bulge mass to satellite mass (val- 
ues 1, 3, 6), and describes the orbit (D for direct, R 
for retrograde). Codes are given in Table ||. On one of 
the merger remnants, we ran a second merger, which we 
denote H3R3R (see § || and inset in Figure ||). We also 
attempted a multi- merger experiment involving 10 small 
satellites (§ ||). 

Computations were carried out with a SGI Power 
Challenge machine (6 64-bit R8000 processors) . Evolution 
was computed using the TREECODE of C. H. Heller (sec 
Heller 1991, Heller & Sholsman 1994), kindly made avail- 
able by the author. Heller's code, with SPH turned off, 
uses the cubic tree structure described by Barnes & Hut 
(1986). The algorithm updates the particle positions with 
the leap-frog algorithm, with a variable time step rang- 
ing between 0.01 and 0.05. The gravitational force was 
softened with a spline kernel (Hernquist & Katz 1989) 
with constant softening length e = 0.02. No quadrupole- 
moment corrections were applied. In cubic treecodes, a 
monopole-only calculation increases force errors by about 
50% relative to a calculation including quadrupole terms 
(Hernquist 1987). Hence, for the tolerance parameter used 
(6 = 0.8), our code computes forces within 1.5% of those 
given by direct summation, compared to typical 1% errors 
of the calculation to quadrupole order. All models were 
evolved beyond the full merger. Total energy was con- 
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served to better than 0.1% during the simulation. Two of 
the merger models were repeated using Hernquist's version 
of the TREECODE (Hernquist 1987, 1990a). The final 
surface density profiles were identical to those computed 
with Heller's code. 



3. Results 

In figures | and | we show snapshots of the evolution of the 
models HID and H6D. In both these models, as in models 
1-6, 8 and 9, the core of the satellite reaches the center of 
the primary galaxy and merges with the bulge. The scaling 
of dynamical friction with the satellite mass causes the 
merger times to significantly increase as the satellite mass 
decreases. The merger affects the structure of the main 
galaxy in several ways. Transient spiral patterns, warps 
and non-axisymmetrical structures appear in the disk of 
the main galaxy. We address those in § [|. 

We first focus on the structural parameters of the 
bulge. Figure ||a shows the surface density profile of the 
main galaxy before the merger. Poisson-based error bars 
in the density, not plotted, are a few percent at most owing 
to the large number of particles. We decompose this pro- 
file in two components much in the same way as observers 
decompose a galaxy surface brightness profiles. The disk 
portion of the profile is well fit with an exponential law. 
The Sersic law written as 



J(r) = I e exp [-(0.868n- 0.142) {( £ f /n - l} 



(1) 



is used for the bulge (Sersic 1968; Andredakis et al 
1995; Prieto et al. 1999), with I e , r e and n as free param- 
eters. The best fit is achieved with n — 1, i.e. exponential 
profile. Deviations occur only at radii smaller than the 
core radius, where the surface density profile flattens. 

Figures ||b-j show the final face-on, azymuthally- 
averaged radial surface density profiles of all luminous 
matter for models 1-9. The dotted lines show the two- 
component, exponential plus Sersic-law fits to the surface 
density. Again these are raw fits to the luminous matter, 
ie. no account is made of whether particles originally be- 
longed to the bulge, to the disk or to the satellite when 
performing the decomposition of the profile. Thus, we an- 
alyze the mass distribution in much the same way as an 
observer would model the light distribution of a face-on 
spiral galaxy. Points inside r = 0.15 were excluded from 
the fit, as they are too sensitive to the flattening of the 
initial bulge surface density profile. The fits were made 
using a Levenberg-Marquardt nonlinear fitting algorithm 
to locate the x 2 minimum. All model parameters were al- 
lowed to change. The entire set of fitted parameters and 
the bulge-to-disk mass ratio derived from the fits are tab- 
ulated in Table |. The final n increases with the satellite 
mass, being slightly higher in the retrograde cases (7%) 
than in the prograde cases, n reaches 4 for a satellite as 
massive as the original bulge. In order to compute the 



goddness of the fits we computed the reduced \ 2 f° r fits 
in which n was fixed at values ranging from 0.8 to 5. Fig- 
ure ||a-j shows the reduced x 2 vs n - arrows indicate the 
minimum \ 2 values listed in Table ^. The minimum val- 
ues of x 2 together with the number of data points in the 
fit n = 30, imply a better than 99% confidence that the 
errors are not due to mismatch of the model to the data. 
Figure |^a-j demonstrates that the remnant bulges have 
indeed evolved away from the n = 1 initial exponential 
profile, given that values of n below the best-fit values are 
strongly excluded. The x 2 profiles have well-defined min- 
ima, with the exception of the massive-satellite models 
HID, H1R, KIR. These models admit fits with n > 4. 

The models displayed in figures ||c-h all correspond to 
Hernquist satellites. Models H1R and H3R were repeated 
using King model satellites, as a basic test of the depen- 
dency of the result on the choice of satellite model. These 
are models 8 (KIR) and 9 (K3R) in Table §. Masses and 
half-mass radii for the satellites of models KIR and K3R 
are identical to those of models H1R and H3R, respec- 
tively. The surface density profiles of models KIR and 
K3R (Fig. ||,j) are very similar to those of H1R, H3R, ex- 
cept at the very center. The Sersic shape index n, given in 
each figure, is nearly identical in models KIR and K3R as 
in models H1R and H3R, respectively. This suggests that 
the present results do not depend on a fundamental way 
on the details of the internal structure of the satellite, but 
mostly on the satellite total mass and mean density. 

Figure ^| shows growth vectors in the n-\og(B/D) 
plane. As discussed above, the bulge-to-disk ratios are 
those derived from the Sersic plus exponential fits, i.e. 
they are the ones an observer would derive. Overall, bulges 
grow with the accretion process. As satellite accretion de- 
posits mass onto the central bulge, it drives the increase of 
the n shape parameter of the surface density profile. Plot- 
ted in the same figure is the distribution of n vs B/D for 
the bulges in the BP94 sample, extracted from Figure 5b 
of APB95. The match is encouraging, as the growth vec- 
tors trace an increase of n Sersic index with B/D just as 
observed in the data; and most of the range of n values 
displayed by the data is obtained in the models. Admit- 
edly our models cover the high-mass end of the B/D space 
only. Also, the growth vectors trace a somewhat steeper 
slope than that of the distribution of data points. Indeed, 
accretion appears as being effective in driving bulges to- 
ward high values of n. Two conclusions can be drawn from 
this result. First, collisionless accretion of dense satellites 
onto disk-bulge-halo galaxies causes the bulge surface den- 
sity profile to evolve toward higher-n Sersic profiles. That 
a merger of similar-mass systems leads to an r 1 / 4 profile 
is well known (Gerhard 1981, Barnes 1988). The present 
mergers involve very unequal masses, thus our result is not 
a direct consequence of theirs, and indeed the final (total) 
mass distribution in our models differs substantially from 
the r 1 / 4 law. The second conclusion is that bulges showing 
exponential surface brightness profiles may not have sig- 
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Fig. 3. Radial surface density distribution of luminous matter. Error bars are smaller than the symbol size, (a) 
Initial model, (b) Low-density model 7 after the merger is complete, (c-h) Models 1-6 after the merger is complete. 
Points: model measurements. Dotted lines: Exponential and r 1 /™ fitted components. Solid line: Sum of the two 
fitted components. 



nificantly grown via collisionless accretion of dense satel- 
lites. 

As said in § |^, we let the remnant of model H3D absorb 
a second satellite of the same mass as the first. The inset to 
Figure ^ shows the growth vector of the resulting model, 
H3D3D, together with those of H3D and HID. The second 
merger results in a similar increase in n as the first merger, 
as well as a similar fractional growth of B/D. The figure 



suggests that a subsequent merger with a similar satellite 
(total accreted mass in the three mergers equal to the 
initial bulge mass) would bring the final remnant to a 
n and B/D comparable to those of model HID (single 
satellite with mass equal to the initial bulge mass). This 
suggests that the evolution of bulges in the n-B/D plane 
has little dependence on whether the mass is accreted in 
one event or peacemeal. 
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Fig. 4. Radial surface density distributions of The various luminous components, after the merger. Asterisks: bulge 
material. Diamonds: disk material. Dotted line: Satellite material. Near the x-axis origin, Thick arrow: Half mass 
radius of the bulge matter distribution. Thin arrow: Half mass radius of the satellite matter distribution. 



The satellite needs to reach the galaxy center undis- 
rupted for it to modify the bulge shape parameter. In the 
low density satellite case (model 7, LID), the satellite dis- 
rupted completely before reaching the bulge. The decom- 
position of the surface brightness profile in this case is 
shown in Figure ||b. The best fit for the bulge has n = 1, 
the same value as in the initial profile. The mass distribu- 
tion is however not entirely undisturbed. Comparison with 
the initial surface density profile (Fig. H]a) shows that the 



final extrapolated disk central surface density of model 
LID is slightly lower, and, more significantly, the central 
density of the bulge component is higher than in the ini- 
tial model. This is surprising as the satellite has disrupted 
completely before reaching the center. 

We now analyze the contribution of each mass compo- 
nent to the final surface density profiles. Figure |4| shows 
the surface density distributions for matter originally be- 
longing to bulge, disk and satellite, for all of the merger 
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Fig. 5. Reduced x 2 values for the diffcrents fits. The mininum value of x 2 is indicated by arrows in each plot. 



models. Panel (a) gives the distributions in the models 
before the merger. Near the x-axis of each plot, vertical 
arrows mark the half-mass radius ri/ 2 for the distribu- 
tion of bulge particles (thick arrows) and that of satellite 
particles (thin arrows). 

The distributions of particles initially belonging to the 
bulge and to the satellite develop extended tails in all 
cases, with ri/2 being larger for the more massive satellite 
cases. For the bulge, the tail is a result of the absorption 
by dynamical friction of the orbital energy and angular 
momentum of the satellite, while for the satellite it is a 



result of stripping. Matter from both the bulge and the 
satellite contribute to the final profile of the inner compo- 
nent deviating from the straight, exponential shape. For 
the more massive satellite mergers, the central density of 
bulge material drops significantly, the center being filled 
with satellite material; the central density of bulge parti- 
cles is roughly unchanged for lower mass satellites. 

There is a third contribution to the final shape of the 
inner component deviating from an exponential shape - 
that of the disk material. In all cases, the surface den- 
sity profile of particles originally belonging to the disk 
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significantly deviates from an exponential law, curving up 
inward and showing a significantly higher central surface 
density than in the initial model. Disk material has been 
dragged inward during the merger. When performing raw 
two-component fits such as those shown in Figure |^, mat- 
ter originally belonging to the disk contributes to the in- 
ner component. This occurs even in model LID where the 
satellite has completely disrupted before reaching the cen- 
ter. 

How much damage does the initial bulge take as as 
result of the absorption of orbital energy? Figure ^ gives 
R m , the radius enclosing a given fraction of the mass, 
for the bulge particles in the initial and the final snap- 
shots for models 1-7. In models 1 (HID) and 2 (H1R), the 
initial bulge expands at all radii, reflecting the effects of 
the denser satellite on the relatively fragile bulge. Smaller 
satellites (models 3-7) deposit most of their energy in the 
outer layers. In all cases, the increase in R m traces the 
formation of an extended tail which turns the initial ex- 
ponential profile into an n > 1 Sersic profile. 




-1 

log(B/D) 



4. Internal kinematics of the merger remnants 

The transformation brought forward by the merger modi- 
fies not only the galaxy's radial mass distribution but also 
its internal kinematics. Rotation curves of disk galaxies 
follow systematic patterns as a function of Hubble type, 
hence we may obtain additional diagnostics on the accre- 
tion by studying the rotation curves of the merger rem- 
nants. The rotation curves of the initial model and of mod- 
els 1-9 after the merger are shown in Figure 0. These are 
line-of-sight velocities, measured along a virtual slit placed 
in the plane of the disk, with the disk seen edge-on. Veloc- 
ities for particles initially belonging to each of the galax- 
ian components are shown with different symbols. For the 
post-merger rotation curves, bulge and satellite particles 
have been grouped together for clarity. The individual ro- 
tation curves for bulge and satellite material are shown in 
Figure |[ 

We first analyze the direct models (Models 1, 3, 5, 7). 
The overall rotation curve of the luminous matter (Fig- 
ure 0, solid line) is steeper in the model representing the 
earliest type (model HID, largest final bulge), becoming 
shallower for later types (model H6D). This behavior is 
similar to that observed in disk galaxies (eg. Casertano 
& van Gorkom 1991). The faster rotation of the model 
with the largest bulge reflects not only the contribution 
of the massive bulge to the gravitational potential but 
also the intrinsic fast rotation of the bulge. In the cen- 
tral parts this is higher than that of the disk material, 
whose motion is partially pressure-supported due to heat- 
ing: merger 1 has produced an object that may resemble 
an SO galaxy. It is the material originally belonging to 
the bulge that drives the fast inner rotation of the final 
bulge (Fig. |8^,), while the satellite material rotates much 
slower. The difference is due to the differences in the den- 



Fig. 6. Growth vectors in the n-\og(B/D) plane. Each 
arrow starts at the location of the original model and ends 
at the n and B/D derived from the two-component fit 
to the surface density profile after the merger. Growth 
vectors for models H3D and H3D3D have been displaced 
to the left for clarity. They are displayed in arbitrary units. 
Points: The distribution of n vs. \og(B/D) for bulges, 
from Fig. 5b of APB95. 



sity profile of bulge and satellite. The bulge (King model) 
and satellite (Hernquist model), despite having equal total 
masses, have different central densities; when both merge, 
the satellite high-density core eventually dominates the 
potential, driving bulge material to spin in the satellite's 
wake. This behavior is common to all mergers of galaxies 
with unequal densities (Balcells & Quinn 1990, Balcells & 
Gonzalez 1998). In models involving King satellites, the 
rotation curves of bulge and satellite materials are more 
nearly similar (Fig. ||i,j). Smaller satellites have overall 
lower effect on the total potential. Stripped, high veloc- 
ity material dominates the rotation curve of the material 
originally belonging to the satellite (Fig. |]g,h). 

Retrograde models (nos. 2, 4, 6, 8, 9) show characteris- 
tic retrograde rotation of the bulge and satellite material. 
In the central region, the final bulge spins due to deposi- 
tion of orbital angular momentum transported inward by 
the satellite, and the material belonging to the initial bulge 
acquires rotation (our initial models featured non-rotating 
bulges). Outside the bulge (at roughly r > 0.5), stripped 
satellite particles and particles initially belonging to the 
bulge orbit as test particles in the halo potential. This 
material may correspond to the retrograde moving groups 
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Fig. 7. Line-of-sight velocities of the models, observed from a viewing angle in the plane of the disk, (a) Initial model. 
(b) Low-density model 7 after the merger is complete, (c-h) Models 1-6 after the merger is complete. Stars: disk 
material. Dashed line: Dark matter halo. Solid line: Luminous matter. 



in the Milky Way halo (Norris & Ryan 1989, Majewski 
1992), but would be unobservable in unresolved galaxies 
due to its low surface brightness. 

The counterrotation of the final bulge is common to all 
our retrograde mergers. If the initial bulges had had direct 
rotation, we expect that the final bulge would not counter- 
rotate in the low-mass mergers (models H3R, H6R, K3R, 
Fig. |?]f,h,j), but new models are needed to verify this. For 
the low-mass retrograde mergers, the fast counterrotating 
velocity of the satellite material (Fig. |sja) , can be detected 



with line-profile spectral analysis techniques such as un- 
resolved Gaussian decomposition (Kuijken & Merrifield 
1993). These simulations show that the counterrotation 
found in the bulge of NGC 7331 (Prada et al. 1996) could 
correspond to the remnant of an accreted low-mass satel- 
lite. Bulge counterrotation has also been found in NGC 
2841 (Silchenko et al 1997), but is otherwise uncommon 
in spiral and SO galaxies (Kuijken et al. 1996). 
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5. Effects on the disk 

It is well known that mergers with dense satellites heat up 
the disk (QHF93), and the effect is clearly seen in Figures 
| and | Figure [id shows the disk scale height zjj for all 
models as a function of time, zjj has been measured at the 
galactocentric distance of R — 1 which corresponds to one 
initial disk scale length. The thickness of the disk increases 
with the mass of the satellite galaxy. Satellites with the 
same mass as that of the bulge increase zd by about 2.5 
times more than satellites with 1/3 of the bulge mass. No 
significant differences can be found between prograde an 
retrograde models. Figure [II] gives also the thickness of 
the disk for the low-density satellite model (L1R). The 
increment of zd in model L1R is about 3.4 times less than 
the increment presented by models HID and H1R in which 
the satellite has the same mass but higher density. Our 
thickening results arc comparable to those of QHF93. 

A second important effect of the merger is to increase 
the disk scale length, h increases by 10% for the smaller, 
increasing to 60% for the largest satellites (Table ||) . The 
shallower disk profile is a result of outward transport of 
disk material in the outer parts, combined with inward 
transport to the bulge in the inner parts. Because the 
bulge effective radius does not increase strongly during 
the merger, the ratio r e /h decreases as a result of the 
merger. The final r e /h does not scale with the final B/D 
in any systematic way. We discuss the significance of these 
results in § ^. 

6. Discussion 

Our simulations show that, as a spiral galaxy grows its 
bulge by accretion of small satellites, the bulge surface 
brightness profile quickly evolves from an n = 1 , exponen- 
tial profile, to a profile approaching n — 4. Growth vectors 
on the n vs B/D plane suggest that the dependency of n 
on Hubble type found by APB95 could be the result of 
satellite accretion. 

The result goes beyond being just a manifestation of 
the known evolution of violently relaxing systems toward 
the r 1 / 4 law (eg. van Albada 1982). Indeed, gravitational 
matter in our systems includes disk and halo in addition 
to the bulge. Taken as a whole, none of the final rem- 
nants approach the r 1 / 4 law, especially so the luminous 
components. Our results demonstrate that, only that spe- 
cific subsystem identified with the central bulge, gradually 
evolves toward r 1 / 4 , while the rest of the luminous matter 
keeps its exponential surface density distribution. 

Collisionless processes alone are involved in the trans- 
formation. Satellites that do not entirely disrupt during 
the merger are needed for the process to operate with the 
efficiency shown here. If the satellite disrupts the bulge 
surface brightness profile remains undisturbed. However, 
the transformation does not rely on the satellite having 
an r 1 / 4 to begin with: the evolution of the surface bright- 
ness profile is similar whether the satellite is modeled with 



a Hernquist profile or a King profile, suggesting that, as 
long as the satellite does not disrupt during the merger, 
accretion-driven bulge growth makes the bulge evolve from 
n = 1 toward n — 4 profile shapes. Indeed, the evolution 
of the bulge surface brightness profiles is driven more by 
the puffing up of the bulge material by the absorption of 
orbital energy and angular momentum of the satellite than 
by the deposition of the satellite's high-density cusp in the 
remnant center. We expect that the details of the central 
surface brightness profiles do depend on the shape of the 
central density profile of the satellite (eg. compare the cen- 
tral profiles of H1R to KIR, and H3R to K3R, Fig. |), but 
the effects on n are small. We do not make quantitative 
predictions on the resulting central densities because our 
King models deviate from exponentials within r = 0.16, 
and because softening limits the ability of the models to 
accurately reproduce central densities. 



ft; 




2 4 6 8 10 

Model 



Fig. 9. Radii enclosing a given % of the mass for the dis- 
tribution of particles initially belonging to the bulge. The 
abscissa is the model number from Table ||. Model rep- 
resents the initial bulge set of R m . The horizontal dotted 
line is the radius initially enclosing 50% of the bulge mass. 



Mixing of the initial populations occurs not only due to 
heating but also because the radial redistribution of disk 
material toward the center makes some material initially 
belonging to the disk to be interpreted by the observer 
as belonging to the bulge when fitting the disk with a 
standard exponential profile. This mixing may contribute 
to the extreme color similarity between bulges and inner 
disks (Terndrup et al. 1994, Peletier & Balcells 1996). Mix- 
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Fig. 8. Line-of-sight velocities of the particles originally in the bulge (full line) and in the satellite (diamonts) after 
the merger is complete, from a viewing angle in the plane of the disk. 



ing does not occur at the very center, where the satellite 
core ends up, hence on purely stellar dynamical grounds 
we should expect a color signature there. Detailed inter- 
pretation however is made uncertain due to the presence 
of dust and star forming processes. The color structure 
in the inner few 100 pc of bulges is now known to be far 
from simple, with steep reddening profiles and hidden star 
formation processes revealed by HST/NICM0S-WFPC2 
imaging (Peletier et al. 1999). Gas in the disk, not mod- 
eled in our simulations, is likely to respond more violently 



to the accretion and pile up at the center (Barnes & Hern- 
quist 1996) probably triggering star formation and chem- 
ical enrichment. 

The results are somewhat affected by the choices made 
for the simulation, such as the bulge rotation, the satel- 
lite rotation or the types of merger orbit. Rotation in the 
initial bulge, not included in our models, would probably 
make the prograde and the retrograde cases more different 
from each other, broadening the range of evolution vectors 
in Figure [j[ Otherwise, rotation is not likely to signifi- 
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Table 2. Orbital parameters for the merger experiments 



Model 
(1) 


Code 
(2) 


Ms at /Mb 
(3) 


bat 
1 1/2 

(4) 


V r 
(5) 


V e 
(6) 


Inclin (deg) 
(7) 


1 


HID 


1 


0.243 


-0.001 


0.222 


30 


2 


H1R 


1 


0.243 


-0.001 


0.222 


150 


3 


H3D 


0.333 


0.102 


-0.025 


0.227 


30 


4 


H3R 


0.333 


0.102 


-0.025 


0.227 


150 


5 


H6D 


0.167 


0.068 


-0.115 


0.219 


30 


6 


H6R 


0.167 


0.068 


-0.115 


0.219 


150 


7 


LID 


1 


0.347 


-0.011 


0.222 


30 


8 


KIR 


1 


0.243 


-0.001 


0.222 


150 


9 


K3R 


0.333 


0.102 


-0.025 


0.227 


150 



Description of the columns: (1) Model number. (2) Model code. (3) Initial mass ratio between satellite and bulge. (4) Initial half-mass radius of the satellite. (5) and (6) 
Radial and tangential velocity components of the relative orbit. (7) Initial angle between the orbital angular momentum and the disk spin. 
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Fig. 10. Time evolution of the disk scale height zd for the 
models. This scale length was measured at R=l (see test 
for details). 



cantly affect the evolution of n. Differences may occur if 
the bulge, or the satellite, had a central black hole. Mas- 
sive black holes are now virtually confirmed in a number 
of galaxies, and they might be by-products of the initial 
formation of the bulge (Rees 1993). Central black holes 
would probably affect both the final density profile and 
the rotation of remnant. Modeling disk galaxy mergers 
with smaller initial bulges would be interesting, although 
difficult computationally. 

The lack of bulge initial rotation and the higher density 
of the satellite both contribute to making counterrotation 
more prevalent in the models than they are likely to be 
in reality, and it is not a prediction of our models that 



half of all bulges should contain counterrotating compo- 
nents. Quantifying merger- induced bulge counterrotation 
will require a comprehensive set of simulations that in- 
clude a realistic distribution of satellite orbit orientations 
(Zaritsky et al. 1997). 

The accretion of dense satellites is less effective in mak- 
ing the bulge's effective radius r e grow than it is at making 
the disk scale length ho grow (Table || from absorption 
of orbital energy and angular momentum. The decrease 
of r e /fiD with each merger is at odds with the observed 
increase of r e /ho toward early types in real galaxies (Gra- 
ham & Prieto 1999). Two explanations for this mismatch 
are possible. First we note that the dense satellites de- 
posit most of their mass in the center of the bulge, while 
the low density model (LD1) deposits all of its mass out- 
side the bulge. We may envision that for intermediate den- 
sity satellites the mass is deposited throughout the bulge 
driving a stronger increase in the bulge's r e . Second, we 
discuss below that the surviving disk in in fact a thick 
disk. De Grijs & Peletier (1997) measure the scale lengths 
of thick and thin disks in edge-on galaxies, showing that 
hthick » hthin- This may explain that r e /hu decreases 
in our models, where hjj = H t hick, while it increases in 
real galaxies, where hp = h t hin- 

The disks thicken during the process. The exact val- 
ues shown in Fig. 10 are sensitive to the choice of energy 
and inclination of the initial orbits, with higher thicken- 
ing expected of more inclined, and more energetic, or- 
bits. Discussion of the disk vertical heating effects are 
beyond the scope of this paper, and are extensively ad- 
dressed in eg. QHF93. Our thickening results are similar 
to theirs, and may be interpreted as indicating that the 
merger contributes to the growth of a thick disk. Lower- 
ing the satellite mass to reduce the impact on the disk 
helps a bit. Our low- mass satellites impart less damage 
to the disk. But their effect on the surface density of the 
bulge is also limited (n — 1 to n ~ 1.5). It is plausible, 
but not obvious a priori, that a sequence of such mergers 
may drive n to keep growing. Lowering the satellite mass 
even further is of limited use, as dynamical friction be- 
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Description of the 
Bulge- to- disk final 



Model 


Code 




Bulge 




Disk 


B/D 


(1) 


(2) 


log(Ie) 
(3) 


r e 
(4) 


n 
(5) 


log(Io) 
(6) 


h 

(7) 


(8) 


Initial 




0.04±0.05 


0.21±0.02 


1.05±0.13 


-0.80± 0.03 


1.0±0.03 


0.51 


1 


HID 


-0.05±0.01 


0.23±0.01 


4.03±0.11 


-1.23± 0.02 


1.63±0.08 


1.26 


2 


H1R 


-0.21±0.06 


0.23±0.01 


4.30±0.23 


-1.03± 0.04 


1.20±0.07 


0.92 


3 


H3D 


-0.20±0.03 


0.24±0.05 


2.35±0.15 


-1.12± 0.05 


1.45±0.05 


0.72 


1 


H3R 


-0.18±0.04 


0.22±0.03 


2.51±0.09 


-0.90± 0.03 


1.11±0.06 


0.85 


5 


H6D 


-0.12±0.02 


0.21±0.01 


1.53±0.12 


-0.96± 0.01 


1.19±0.03 


0.53 


6 


H6R 


-0.18±0.03 


0.20±0.02 


1.65±0.16 


-0.92± 0.06 


1.13±0.04 


0.43 


7 


LID 


-0.01±0.01 


0.21±0.02 


1.03±0.11 


-0.89± 0.07 


1.12±0.03 


0.51 


8 


KIR 


-0.05±0.02 


0.23±0.02 


4.01±0.13 


-0.90± 0.03 


1.05±0.02 


1.24 


9 


K3R 


-0.10±0.03 


0.23±0.04 


2.53±0.17 


-0.90± 0.02 


1.06±0.03 


0.88 


column*: (3) 


Bulge effect 
derived fror 


ivc surface dent 
a the Sersic plu: 


iity, (4) Bulge effective radius, (5) Bu 
; exponential double fit. 


Igc profile index, (6) 


Disk central inten 


sity, (7) Di 



scale length, (8) 



comes inefficient. Satellites with mass 0.1 times the bulge 
mass have merger times approaching one Hubble time. We 
ran a multi-merger case in which 10 satellites with mass 
0.1 times the bulge mass evolved around the disk galaxy. 
Weak dynamical friction and satellite-satellite interactions 
prevented the merger over one Hubble time. 

Therefore, disk thickening is probably unavoidable, 
and, because we start out with an already formed disk, the 
present scenario of bulge growth is variation on the " bulge 
before the disk" theme. Indeed, for such galaxy to match 
present day disk galaxies, a thin disk must be rebuilt out 
of the gas remaining after the last merger. The formation 
of the thick disk is itself not a problem for this scenario, 
as thick disks are present in early type disk galaxies and 
absent in late- types (de Grijs & Peletier 1997). 

Merger-driven bulge formation prior to the disk (e.g. 
Kauffmann et al. 1994) is at odds with the fact that Sc 
and later-type spirals have exponential bulges: our experi- 
ments indicate that exponential bulges are extremely frag- 
ile against accretion and merging. The seeds of the central 
bulges of spiral galaxies must be given by events unrelated 
to accretion of smaller galaxies. 

The efficiency of the increase of n (Fig. ^|) places limits 
on the fractional growth of n = 1 — 2 bulges by collisionless 
accretion of dense satellites. Our results suggest that col- 
lisionless accretion of dense satellites bear the signature 
of a non-exponential bulge surface density profile. Such 
signature lends itself to tests requiring imaging data only, 
hence can be be used to constrain the merger history of 
galaxies at high cosmological distances. 

7. Summary 

The collisionless accretion of dense satellites onto disk 
galaxies drives the growth of the bulge and an increase 
in the n index of the Sersic fit to the bulge surface density 
profile. The mass of the bulge and n grow proportional 
to the satellite mass. A single merger with a satellite as 
massive as the bulge forms an r 1 / 4 bulge. The range of 



n values obtained in the models matches that observed 
in bulges by APB95, which points at accretion as a sim- 
ple way of setting up the APB95 relation. These results 
support the idea that bulges of late type spirals which 
show exponential surface brightness profiles can evolve to 
bulges of early type galaxies with n > 1 by mergers with 
satellites galaxies. 

We predict a fair amount of population mixing by ex- 
pansion of the bulge material and piling up of disk mate- 
rial to the center. 

The rotation curves of the merger remnants are steeper 
for more massive bulges, in accordance with observations. 
Retrograde mergers may lead to counterrotating bulges. 

The models predict that thick disks form as by- 
products of the evolution toward larger B/D, in accor- 
dance with the presence of thick disks in SO and Sa disk 
galaxies. Matching to present day galaxies requires the 
rebuilding of a thin disk out of remaining gas. 

The efficient transformation of exponential bulges by 
accretion suggests that collisionless mergers at high-z are 
not responsible for the structure of bulges of late-type spi- 
rals. 
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